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We investigate how short and long electron trajectory contributions to high harmonic emission and their in- 
terferences give access to intra-molecular dynamics. In the case of unaligned molecules, we show experimental 
evidences that the long trajectory signature is more dependent upon the molecule than the short one, pro- 
viding a high sensitivity to cation nuclear dynamics within 100's of as to few fs. Using theoretical approaches 
based on Strong Field Approximation and Time Dependent Schodinger Equation, we examine how quantum 
path interferences encode electronic motion whilst molecules are aligned. We show that the interferences are 
dependent on channels superposition and upon which ionisation channel is involved. In particular, quantum 
path interferences encodes electronic migration signature while coupling between channels is allowed by the 
laser field. Hence, molecular quantum path interferences is a promising method for Attosecond Spectroscopy, 
allowing the resolution of ultra-fast charge migration in molecules after ionisation in a self-referenced manner. 



I. INTRODUCTION 

In the intuitive description of the high-harmonic gener- 
atio n (H HG) mechanism using a linearly polarised laser 
fieloP^ an initially bound wave packet is first freed in the 
continuum near the peak of the laser electric field (step 
1: ionisation). During this step the system consists of 
a superposition of a continuum electronic wave packet 
(EWP) subject to the laser field and an EWP remaining 
in its ground state. The continuum EWP follows a tra- 
jectory where it is accelerated away and returned to the 
core driving by the oscillating laser field (step 2: propa- 
gation). At the returning time, this EWP can then in- 
terfere with the part that remained bound in the system 
and release the extra energy gained during the trajec- 
tory into a coherent harmonic photons emission (step 3: 
recombination). The harmonic photons emitted cover a 
large range of frequencies from UV to XUV and the cor- 
responding spectra presents a specific extended comb of 
harmonic frequencies with a plateau followed by a cut- 
off region. This large spectral range allows experimen- 
tally the production of an isolated attosecond pulse^E! 
or a train of XUV attosecond pulses^. Classically, the 
instant at which the EWP is freed in the continuum, re- 
ferred to as ionisation time, determines the trajectory the 
continuum EWP will follow durin g the propagation step. 

It is now well-establishe d 2 ! 7 ! 8 ! that in the plateau re- 
gion multiple families of trajectories can contribute to 
the same harmonic emission whilst in the cut-off (highest 
photon energies) one trajectory family contributes only, 
referred to as the 'cut-off trajectory'. To give a first es- 
timation of these trajectories one can simply resolve the 
classical Newton equation of motion of a charge particle 
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(electron) in a laser field over one single laser cycle as de- 
scribed in figure 1. The three type of trajectories are pre- 
sented in l.a with a colour map representing the kinetic 
energy the electron acquired among a specific trajectory. 
This kinetic energy is then plotted as a function of ion- 
isation and recombination time in l.b. The maximum 
kinetic energy, which corresponds to the highest photon 
energy emitted, results from the cut-off trajectory. From 
this trajectory one can define within an optical cycle two 
families of trajectories, referred to as the 'short' (with 
duration smaller that the cut-off one) and the the 'long' 
(with duration longer that the cut-off one). The short 
trajectories are characterised by a kinetic energy that in- 
creases with the excursion time while the long trajectory 
has a kinetic energy that decreases as a function of the 
excursion time in the continuum. In both case this results 
in an attosecond chirp (positive for the short trajectory 
and negative for the long) that provides a time-frequency 
mapping. Therefore dynamical effects of the core (time 
domain) during HHG are encoded in the high harmonics 
spectra (frequency domain). For kinetic energy below the 
cut-off energy (maximum) , there exists therefore a short 
and a long trajectory that both lead to the emission of 
the same photon energy. These are the two trajectories 
that can interfere leading to the method we propose here 
for studying Attosecond spectroscopy. 

In the quantum picture and under the strong field ap- 
proximation these trajectories correspond to Feynmans 
paths where the harmonic dipole moment of a given har- 
monic order q can be written as a Fourier sum over dif- 
ferent quantum paths j contributing to the XUV emis- 
sion. Hence the quantum paths are a generalization of 
the semi-classical electron trajectories described in the 
three steps model. The time dependent dipole induced 
is given as follow: 
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FIG. 1. Classical calculations: Left hand l.a: trajectories as 
a function of the ionisation and returning time. Right hand 
l.b: corresponding kinetic energy acquired by the electron 
freed in the laser field. Cut-off trajectory (red) corresponds 
to the highest kinetic energy leading to the so called cut-off 
law (1.3/p + 3.2f/ p ). At a given kinetic energy below the cut- 
off threshold, two trajectories referred to as 'short' and 'long', 
contribute to the same kinetic energy 
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where q is the harmonic order, j define the quantum 
path, dq corresponds to the dipole amplitude over a 
quantum-path j and $^ is the respective dipole phase: 



*i = W + f^2* (2) 



where I p is the biding energy of the ground state (first 
ionisation potential level), T 3 q is the excursion time also 
known as the duration of a quantum-path j correspond- 
ing to the difference between ionisation times and re- 
combination times, p is the canonical momenta and A 
the vector potential. Because the ionisation time win- 
dow differs from the recombination time window, high 
harmonics are not emitted at the same time^. This lack 
of synchronisation, the well known attosecond chirp also 
present in the classical description, is actually an advan- 
tage for Attosecond spectroscopy as this way each har- 
monic order maps a given time in the laser field. This 
time scale depends on the laser field wavelength Ao and 
the laser peak intensity Iq , as the cut-off position scales 
with Aq and Iq, but also on the trajectories followed. The 
mapping law in the plateau for each trajectories is given 
by the short and long trajectories excursion times (du- 
ration of each path). To estimate this law we performed 
quantum calculations of the harmonic dipole (under the 
Strong Field Approximation) and we derived an analytic 
expression of the EWP excursion time for the short and 
the long trajectories as a function of the harmonic or- 
der q performing a fit (that for the short trajectory is 
consistent with the expression derived in M. Lcin irP): 



T sho rt = ^ hf - x arcos{l ^t^)) C3 + — (- x arcos(l ^r^))^] (3) 
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where, u> is the fundamental laser frequency, I p the 
ionisation potential, U p the quiver energy, q the harmonic 
order and C n the specify constants in Tab I. The results 
Tshort and Ti on g are given in optical cycle units. 

Under our experimental conditions (800 nm and I = 
0.5 — 2 x 10 14 VF.cto~ 2 ), we present the mapping law 
in figure 2. As one can see depending on the laser 
peak intensity Iq (colorbar) the cutoff extends toward 
higher harmonic order as expected. For the maximum 
extension (we employed a maximum peak intensity of 



2 x 10 14: W/cm 2 ) the short trajectory maps a time scale 
from ~ 200 asec to the position of the excursion time 
of the cutoff 1.7 fsec whereas the long trajectory maps 
a time scale from the cutoff excursion time (1.7 fsec) 
to ~ 2.7 fsec. Thus, if any intra-cation dynamics oc- 
curs within this time window, the total dipole moment is 
modified dictating that the HHG radiation enc odes key 
signatures of this motion just after ionisatiorPEJ] We will 
indeed examine the case of fast nuclear motion, channel 
superposition signature and channels dynamical popula- 
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TABLE I. Constants for short and long excursion times cal- 
culation 

C\ C2 Cz C4 C5 
short 2 2.1 1.18 0.21 1.06 
long 1.46 2.03 1.53 0.16 1.06 

tions transfer due to laser coupling that happen during 
this time scale. 
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FIG. 2. Excursion time calculated for each trajectories within 
laser peak intensity of 7o = 0.5 — 2 x 10 14 W.cm~ 2 as employed 
in the experiment and shown in the colour bar. The mapping 
law for both trajectories allows us to follow dynamical process 
from 200 asec to 1.7 fsec for the short trajectory and from 1.7 
fsec to 2.7 fsec for the long trajectory 

To access the encoded information, experimental meth- 
ods based on HHG have been employed. The PACER 
technique (Probing Attosecond motion by Chirp En- 
coded Recollision) , that accesses the cation dynamics via 
the time-frequency mapping allowed by the attosecond 
chirp, was first implemented for retrieving nuclear dy- 
namics by comparing the HHG spectrum between dcuter- 
ated ( refere nce) and protenated molecules of the same 
species^Eil The tomography technique^^ permits to 
access information on which ground states (molecular or- 
bitals) is involved in the HHG process and for this a 
normalisation to an atomic partner response (reference) 
is required. Recent works^ based on two HHG source 
interferometry, one source in unaligned molecules (ref- 
erence) and the other one in aligned molecules, shows 
that the HHG signal contains crucial information on the 
electronic orbitals involved in the process and their fast 
dynamics. The ideal tool would be to access these dy- 
namics without recourse to an auxiliary reference to ex- 
tend the possibility of such pump-probe scheme to reveal 
in a more versatile way the dynamical information. 

One route to this self-referenced Attosecond Spec- 
troscopy is to take advantage of the two quantum-paths 
(short and long) and especially of their interferences (self- 
referenced technique referred to as QPI for quantum path 



interferences). Their corresponding ionisation and re- 
combination times being specific to a quantum-path, the 
high harmonics emission time are separated by 100's of 
asec and map the full time scale needed to study dynam- 
ical processes. Note that if in the presented results we 
use 800 nm laser field and make a proof of principle of 
QPI for Attospectroscopy, one way to extend even more 
the time scale of study (see eq.3 and 4) will be to take ad- 
vantages of NIR sources that provide wavelengths (1200- 
2000nm) higher than 800 nm. Hence the molecular QPI 
technique can be extended to more complex molecules 
where dynamics due to ionisation and subsequent struc- 
tural changes will occur beyond 2-3 fsec. 

Preliminary studies on interferences of electron tra- 
jectories involved in the HHG process (know as atomic 
Quantum-Path Interferences QPI) has demonstrated the 
possi bility to access experimentally dipole information in 
atoms^E^. As shown on figure 3, one can see how the 
harmonic yield is modulated corresponding to subsequent 
constructive and destructive interferences of the quantum 
paths with respect to the laser peak intensity. The exper- 
imental data shown in red are consistent with the SEA 
calculations (blue) whereas the TDSE calculation are not 
showing clear modulations with a high enough contrast 
due to the non-selection of the two first quantum-paths 
(short and long) as allowed within the SEA model. The 
periodicity of the modulation close to 0.3 x 10 14 W.cto _1 
is a good indication that the interferences pattern is due 
to short and long trajectories interferences and not higher 
order trajectories (multiple return). 
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FIG. 3. Quantum path interference calculated with TDSE 
(green), SEA (blue) and measured experimentally in Argon 
jet target (red) for harmonic 21 st . 

These studies demonstrated that experimental condi- 
tions can be found where the short and long trajectories 
may interfere, creating an 'interferometer' with a reso- 
lution of 10's asec. The interferometer arms are experi- 
mentally controllable in phase and amplitude. The phase 
of each arm is proportional to the ponderomotive energy 
U p (i. e. depends upon the laser peak intensity or the 
laser wavelength U p ~ I x A 2 ) and the excursion time of 
the electronic wave packet in the continuum specific to 
an harmonic order and a trajectory. Therefore one way 
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to control the relative phase of the interferometer arms is 
by controlling the laser intensity. The relative amplitude 
of the short and long trajectories that contribute to the 
macroscopic HHG signal can be controlled by: adapt- 
ing phase matching conditions, controlling the focusing 
condition of the laser in the gas jetP or using a second 
harmonic of the laser fields. 

Using the phase matching control procedure, both tra- 
jectories can be phase-matched together with different 
divergences which provides a direct spatial mapping of 
the trajectory contributions in the far field. The diver- 
gence of the short trajectory contribution is smaller than 
the divergence of the long trajectory contribution. So the 
short trajectory appears mainly on-axis whilst the long 
trajectory appears off-axis. In the area where both con- 
tribution overlap their interference leads as spatial rings 
forming in the far-field intensity distribution. A far-field 
spatial selection was used to allow the balanc e of t he two 
contributions to maximise the QPI contraslPEH. This 
interferometric measurement offers a high sensitivity to 
any change in the cation which make it suitable for the 
investigation of ultra- fast intra-molecular dynamics. Im- 
portantly QPI is a self-referencing method since it uses 
the observation of the spatial and the spectral distribu- 
tion of the harmonic emission in a single experiment to 
extract the information. 

In this article we show for the first time, an extension 
of the QPI to the case of molecules in order to retrieve 
electronic and nuclear dynamical cation information. 



II. EXPERIMENTAL INVESTIGATION OF NUCLEAR 
MOTION 

All experimental results presented here were performed 
using the same configuration: we used a Thsapphire CPA 
laser system (Red Dragon, KML Inc.) delivering 30 fsec, 
800 nm pulses at a repetition rate of 1 kHz with a max- 
imum pulse energy of 8.5 mJ of which we used 3 mJ in 
the full beam before aperturing for our experiment. The 
laser beam entered a vacuum chamber and was focused 
by a 40 cm radius of curvature (ROC) spherical mirror 
into a pulsed jet of 500/xm diameter aperture with 1 bar 
backing pressure. A flat field spectrometer, composed of 
a grazing incidence grating (87°, 1200 gr/mm, ROC=200 
mm) was used to spectrally disperse the HHG radiation 
in one direction whilst preserving the divergence in the 
perpendicular direction. A multi-channel plate (40 mm, 
Photonis Inc.) was placed in the grating focal plan fol- 
lowed by a backside-illuminated CCD, to allow us to mea- 
sure the spatially and spectrally resolved harmonic emis- 
sion. The acceptance angle of our detector, 20 mrad, 
was sufficient to resolve spatially the short and the long 
trajectories. The gas jet was positioned ~ 2 mm before 
the laser focus to phase-match both trajectories where 
the laser beam cross-section diameter was 75/im. An 
aperture at the entrance of the HHG chamber was used 
to optimise the signal from the distinct quantum paths 



and the QPI, in this case the laser beam was apertured 
down to 4 mm so that the maximum energy throughput 
is ~ 0.6toJ for all the experiment. A variable attenua- 
tor (wave plate/polariser) was placed in the generating 
beam to control the laser peak intensity at the position 
of the jet and thus control the relative phase of the two 
trajectories. The acquisition of each harmonic spectrum 
was performed over 10 000 laser shots. 

Randomly aligned molecules were investigated to 
demonstrate a first experimental evidence of the differ- 
ent dependence of the two quantum path contributions 
upon the fast molecular nuclear motion. Whilst the ex- 
periments are readily extended to aligned molecules, the 
purpose of the present study was to show a first proof 
of principle in the capability of the long and short tra- 
jectories to encode dynamical effects in the HHG signal 
recorded. The nuclear dynamics are expected to change 
the HHG emission amplitude through the e voluti on in 
time of the nuclear autocorrelation functio n 19 * 20 . We 
study respectively the case of C0 2 and N 2 that are ex- 
pected to have a slow evolution of the autocorrelation 
function^, the case of 2 (expected to have a interme- 
diate evolution of the autocorrelation function) and the 
case D 2 and H 2 (expected to have a fast evolution of the 
autocorrelation function but different from each other). 
The accessible time window of our experiment is given 
by the difference of excursion time of the electronic wave 
packet in the continuum between the short and long tra- 
jectories. 

In Fig. 4 HHG is shown in the far field for Argon, C0 2 , 
N 2 , 2 , H 2 and D 2 for an intensity at the position of 
the jet of 2 x 1Q 14 W/ cm 2 . For the atom and most of the 
molecules, each harmonic in the plateau presents a simi- 
lar shape: a central peak surrounded by an annular distri- 
bution. These HHG spatial distributions correspond to 
the signature of short trajectory mainly phase-matched 
on-axis (central spot) and long trajectory mainly phase- 
matched off-axis (ring). The inner spatial ring observed 
around the on-axis region, well resolved in the case of 
N 2 and C0 2 , corresponds to the region of spatial inter- 
ferences of the two trajectories. All measurements were 
made maintaining the same geometry, the same backing 
pressure and the same jet position, allowing comparison 
between all the samples. Figure 4 shows a clear differ- 
ence of the recorded signal from one molecular target to 
the other. The HHG signal of N 2 and C0 2 implies that 
if the nuclear motion is slow for these molecules, both 
trajectories 'see' the same cation at ionisation with mini- 
mal changes by the instant of recombination. In the case 
of D 2 and H 2 the results obtained are highly interest- 
ing. Actually being deuterated and protenated partners, 
these two molecules illustrate the most the signature of 
the nuclear motion differences encoded into the spatially 
and spectrally resolved short and long trajectory contri- 
butions. As one can see the HHG signal is highly af- 
fected, indeed the long trajectory contribution is clearly 
diminished, being almost totally suppressed in the H 2 
case while still present in D 2 . This is consistent with the 
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FIG. 4. Far Field spatial distribution of short (on axis) and 
long (off-axis) trajectories generated in Argon (upper) and in 
randomly aligned molecules (N2,C02,02,D2,H2). We use for 
the HHG a 30 fs IR laser pulse at 2 x lQ^W/cm 2 at the jet 
position located at ~ 2 mm before the laser focus and the 
backing pressure is kept at lbar. The horizontal axis corre- 
sponds to the spectral domain, the vertical axis corresponds 
to the harmonic beam divergence and the colour map is the 
harmonics signal in logarithmic scale. The elliptical 'shadow' 
that cut slightly the 17 th order on the left hand side is the 
edge of the MCP detector. 



idea that the intensity of harmonic radiation is depending 
upon the nuclear motion and is modulated by the square 
modulus of the nuclear autocorrelation function given by 
I < ^n(t)\ip c (t) > | 2 , where ip n (respectively tp c ) is the 
wave packet evolution of the neutral (respectively of the 
cation). Due to the presence of the short and long trajec- 
tories, for a given harmonic order, it is then possible to 
follow this nuclear dynamics over broad time scale, here 
200 asec to 2.7 fsec. The autocorrelation function can 
be calculated for various molecules as nicely shown irP4 
For fast nuclear dynamic as for D2, if 2 the autocorrela- 
tion function decrease rapidly over time compare to O2 
and CO2 and iV 2 having a autocorrelation function evolv- 
ing slowly (autocorrelation function almost constant over 
this time scale). Hence our observation in D2 and H2 is 
a clear experimental demonstration that the faster the 
molecule expand the more affected the HHG signal from 
the long trajectory contribution. The case of O2 appears 
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FIG. 5. Far Field spectrally resolved QPI of 19th for the fol- 
lowing unaligned molecules: CO2, N%, O2, ■ In the case ofI?2 
and H2 we can observe the disappearance of the modulation 
which shows that the long trajectory is not contributing the 
QPI pattern. 



to be intermediate to N 2 ; CO2 and D 2 ; H 2 - Even though 
these molecules do not have the same structure and so 
won't show the same harmonic yield, it illustrates plainly 
the fact that the relative strength of short and long tra- 
jectories is sensitive to nuclear motion of the cation with 
attosecond-femtosecond time scale, especially in H 2 and 
D2 where the comparison is straight. 

To investigate further the role of cation nuclear dy- 
namics we performed an intensity scan of the HHG spa- 
tial distribution. The intensity was varied from 0.5 to 
1.5 x 10 14 W/cm 2 using an attenuator. While changing 
the intensity, the annular distribution evolved from the 
on-axis peak to the annular off-axis distribution. To ob- 
tain spectrally resolved Molecular QPI, we integrate the 
signal over a spatial region located from 10 to 12mrad 
from the central axis (where interference between the 
short and long contributions is the strongest), for all the 
molecules and record the harmonic signal at the harmonic 
frequency as a function of the laser peak intensity. The 
signal obtained presented in figure 5 for all the molecular 
samples, shows a fringe contrast that differs between the 
molecules due to the fact that the strength of the long tra- 
jectory is dependent on the molecules. However the peri- 
odicity seems to be conserved at ~ 0.3— 0.4x 1 14 iy / cm 2 , 
consistent with previous observations in atoms^ 1 ^. This 
confirms that SFA calculation, valid for atomic case, can 
be also use to investigate how the QPI are evolving while 
generating harmonic in unaligned or aligned molecules. 

Thus the difference in the long trajectory amplitude 
observed between each molecular sample can only be at- 
tributed to nuclear motion in the cation between ioni- 
sation and recombination. The Molecular QPI modula- 
tion periodicity appears to be unchanged in all molecules 
studied here indicating the robustness of the continuum 
electron behaviour and setting the conditions for the de- 
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tection of different cation evolution. The observation of 
short and long trajectories in a broad range of molecules 
(Fig. 4) and the confirmation of QPI signatures (Fig. 4 and 
5) suggests the use of molecular QPI to quantify cation 
dynamics. We now consider theoretically how electronic 
dynamics should be accessible from such Molecular QPI 
measurement. 



III. THEORETICAL INVESTIGATION OF MOLECULAR 
QPI AND IONISATION CHANNEL DYNAMICS 

We know from previous studieip!122l that in N 2 the par- 
allel alignment favours the contribution to HHG from 
the HOMO and HOMO-2 orbitals (large dipole due to 
orbital symmetry) whereas the contribution of HOMO-1 
is cancelled (no dipole: nether ionisation nor recombi- 
nation is possible due to orbital symmetry). In perpen- 
dicular alignment the HOMO-1 is favoured (large dipole 
due to orbital symmetry) and HOMO and HOMO-2 are 
diminished (small dipole: ionisation and recombination 
are smaller than in parallel alignment) but not entirely 
cancelled. For a given laser intensity, the HOMO domi- 
nates up to the position of its cutoff and the contribution 
of HOMO-1 is then revealed beyond the HOMO cutoff 
where it dominates the signal over HOMO-2 contribu- 
tion. However in the case of N 2 the difference of ionisa- 
tion potential for HOMO and HOMO-1 is only one pho- 
ton energy at 800 nm, therefore tracking a signature of 
this orbital involved in the extended cut-off while rotat- 
ing the molecule is challenging. Additionally this separa- 
tion of the ionisation channels by just one photon energy 
(800 nm) may allowed a transfer of electronic population 
by laser coupling changing all the orbitals contribution 
weight upon few 100's of asec. 

It is needed to find methods not only based on cut-off 
extension of the harmonic spectra to explore this intra- 
cation dynamics. Proving that the Molecular QPI are 
sensitive to nuclear dynamics of the cation, we will in 
this section exploring the capability of Molecular QPI for 
probing multiple orbitals and molecular cation electronic 
dynamics. We first investigated theoretically the Molec- 
ular QPI signal dependence upon different scenarios for 
N 2 like and model H 2 like molecules. 



A. Atomic model for QPI probing ionisation from multiple 
channels 

Our first approach is very simple and consists in un- 
derstanding how the QPI measurement may be sensitive 
to the ionisation channels. For this we have calculated 
the QPI pattern for a Hydrogen-like atom but with the 
ionisation level of N 2 orbitals: HOMO (I p = 15.58ey) 
named channel 1, HOMO-1 (I p = 16.936^) named chan- 
nel 2 and HOMO-2 (I p = 18.78eV) named channel 3. For 
this we use the Strong Field Approximation and calculate 
the dipole over a multi-cycle laser pulse close to experi- 
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FIG. 6. Atom-like QPI pattern of 19 th harmonic calculated 
for a Hydrogen-like molecule with Is symmetry and N2 ioni- 
sation level of its HOMO, HOMO-1 and HOMO-2 

mental conditions (2 x 10 14 W/cm 2 ; Gaussian envelop; 30 
fs FWHM). The strongest contributions being the short 
and the long trajectories at the high harmonic frequency. 
The results is shown in figure 6 where clear QPI mod- 
ulations are resolved for the 19 harmonic H19 and for 
which the periodicity is found to be ss 0.3 x 10 14 W/cm 2 . 
Similar modulations are found for other harmonics or- 
der. As one can see the QPI pattern presents a clear 
shift for the different channels used, this shift increasing 
toward lower intensity with the I p level. If we consider 
that the QPI pattern, showing subsequent maxima that 
correspond to constructive interferences between short 
and long trajectories, follows a modulated function of 27r 
periodicity hence the channel 2 is shifted by §[27r] with 
respect to channel 1 and 7r[27r] for channel 3 with respect 
to channel 1. This is due to the term I° rbltal T^ in the 
phase acquired by the electronic wave packet upon a spe- 
cific trajectory (j) , tunnelling out and recombining to a 
given orbital (single channel) 2 . Increasing the I p also 
induces a decrease of the recombination and ionisation 
probability therefore as observed the QPI yield is lower 
for channel 3. In this very simple approach no orbital 
symmetry is taken into account nor laser coupling, rotat- 
ing the molecule would correspond to an averaging over 
the QPI patterns of channel 1 (HOMO) and channel 2 
(HOMO-1). Therefore one of the signatures of multiple 
orbitals involved in the HHG process onto the QPI signal 
would be mainly a shift in the intensity position of the 
maxima and minima. This atomic-like' approach shows 
even in a very simple manner that if many orbitals are 
involved their signature should be encoded in the Molec- 
ular QPI. 



B. Molecular model for QPI probing dynamical transfer of 
population between ionisation channels 

During the HHG process, QPI should have the 
potential to separate different scenarios: In aligned 
N 2 molecules, single channel starting from HOMO or 
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HOMO-2 orbitals, and their superposition may con- 
tribute to the HHG. In perpendicularly aligned N2 
molecules, the contribution of HOMO-1 arise and the 
contribution of HOMO-2 can be considered negligible 
(less that 5% of the population is in this state^). In 
this case an additional scenario may be considered where 
a transfer of population between HOMO and HOMO-1 
induced by the laser field is allowed. In a one colour laser 
field at 800 nm, this transfer of population is expected 
to be more favourable in N 2 than in CO2 due to the 
energy gap between these channels (one photon energy 
for N2; three photon energy for C0 2 ). Therefore in a 
first quasi-analytical approach we present a study of the 
QPI in N2 molecule with respect to its alignment to the 
generating laser field. To go even deeper in the under- 
standing of such signature we present in a second study 
an ab-initio calculation of a model H2 molecule subject 
to an additional resonant field that induces this transfer 
of population and allowed to study more accurately its 
signature in QPI pattern of each separated channels. 



1. Semi-analytical model for High-Harmonic Generation in 
molecular Nitrogen 



As beautifully demonstrated by Smirnova et 
nitrogen cation (at a fixed internuclear distance taken at 
equilibrium) has two close ionic states, ionisation from 
HOMO (called X) and ionisation from HOMO-1 (called 
A) that can be pictured with Hartree-Fock method as: 



1 >e 



X 2 Zj : (la g ) 2 {la u ) 2 {2a g ) 2 (2a u ) 2 (lir u ) 4 {3a g ) 



I 2 >= A 2 n+ : (la g ) 2 {la u ) 2 (2a g ) 2 (2a u ) 2 (l7r u f(Sa g ) 

(6) 

In the laser field, this two level system (| 1 > and 
I 2 > with E\ and E2 energies) is split in two quasi-static 
states (i.e. polarised states) | 1 > and | 2 > with ei and 
e 2 energies. In the field free basis and at the ionisation 
time ij, these two states can be written as: 



1 > 

2 > 



1 (Afe) -V(U)\ (\ 1 > 

y/&(U) + V*{U) Viti) A(*0 J \\ 2 > 



(7) 

Where A(t t ) = E 2 -e 1 (t i ) = -(E 1 -e 2 (U)) > 0. At t = U 
the molecule is in a superposition of these two quasi-static 
states: ^ = | 1 > +d 2 {ti) | 2 >. By solving the 

Time Dependent Schrodinger Equation (TDSE) we can 
obtain the evolution of the two quasi-static level system 
in the field free basis. Taking the molecular state after 
ionisation, we can calculate the dipole along each trajec- 
tories (traj = short or long) and for each channel (c = 1 
or 2) using SFA calculations. Using the saddle point ap- 
proximation as in 2 ^, the Fourier transform of the dipole 
can be written as : 



d(fi) 



Aioniti) Ap r0 p(ti , tr)A rec (ti, tr) (8) 



(5) 



where the ionisation time ti and the recombination time 
t r are dependent on the frequency f2 corresponding to 
harmonic energy. Hence: 



i,traj ? tr,traj 



■a l tA[ 
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(9) 



where, Ai on {ti) is the ionisation amplitude at time of ion- 



isation ti, we use the ADK ionisation rate exp - 



3 E L (U) 

The propagation amplitude between the ionisation is 
given by: 



±prop 



(tj, £ r ) 



1 



-i f*T dt'£{t')-il p (ti-t r ) 



x e 



iQ,t T 



(U - t t ) 3/2 

(10) 

And A rec (ti,t r ) =<g\p \ ^ion^k(ti,t r ) > the recom- 
bination amplitude at time of recollision t r between the 
ground state | g > and the exited system ( ^k(ti,t r ) 
being the continuum state and ^i on the ionic state); 



-if*!-dt'at') 



1 ( 2 ^p) 3/2 , with = \p(t'f 



is corresponding to the Volkov propagator 



To take into account the difference of symmetry be- 
tween HOMO (even function along the electron coor- 
dinate; x-axis) and HOMO-1 (odd function along the 
electron coordinate) and the symmetry of the transition 
(without making any assumption on the dipole) we use 
the plane wave approximation of the continuum state 
exp(ifca;) and if we define ^ H for HOMO and vj/^- 1 for 
HOMO-1 as the spatial wave-functions of the ionic states 
then: 



Al\> = / dxV 



s c 

> H (x,y,z)(-i — )exp(ikx) = k.FT x (^ H ) = k / dx^H (x, y, z)cos(kx) = -1 x /i(fc(i;)) (11) 



A l2 > = dx^ H ~ 1 (x,y,z)(-i — )cxp(ikx) = k.FT^ 11 - 1 ) = ik \ dxfH - l(x,y, z)sin(kx) = i x f 2 (k(U)) (12) 



with fi(k) being an even function and /2(A) being an 
odd function. Both functions are set positive for k < 
(the sign of k changing together with the direction of 
laser field). Therefore to represent the even and odd 
characteristics of /1 and j^2^ we have taken: fi(k(ti)) = 
1 if the laser field E(t) > 0, else fi(k(U)) = -1; and 
h(k{U)) = 1. 

For the ionisation amplitude, we need to take into 
account the HOMO and HOMO-1 orbital coordinates 
from which we can determine what is the sign of the 
orbital lobe in the direction of the ionisation to add 
a 7r shift to the ionisation amplitude. The Amplitude 
of I 1 > and I 2 > are both negative if we ionised 
from positive coordinates. If we ionise form negative 
coordinate, the amplitude for | 1 > gets negative and 
the amplitude for | 2 > gets positive. To represent this 

coordinates effect we choose: A^(tj) = — 1 x A^p K (ti) 

■,Af>(U) = lxA% K (U),XE(t)>0 

(ti) = -I xA%> K ( ti ), if E(t)<0 
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FIG. 7. Molecular QPI pattern for 19th harmonic calculated 
for N% as a function of the alignment angle a. The contribu- 
tion of the two channels HOMO and HOMO-1 only differ from 
the single HOMO channel contribution while the molecule is 
perpendicular to the laser 



The transition between the two quasi-static ionic states 
is only permitted in perpendicular configuration of the 
laser field versus the molecular internuclear axis (an- 
gle a) due to molecular orbitals symmetry in N2. In 
other words, the transition dipole being maximum (~ 
0.25a.uf^) for molecules perpendicular to the laser field 
and null for molecules parallel to the field, the par- 
allel component of the coupling is null and the per- 
pendicular component is given by: V(ti,a) = — d.Ei 
= 0.25.Ex(ti) x sin(a)). From this semi-analytic ap- 
proach we calculate the QPI pattern as a function of the 
alignment angle. Results are shown in figure 6 for the 
harmonic 21 st and calculating over one single laser pe- 
riod to reduce the calculation to the short and the long 
trajectories only. 

In the case of a single channel shown in the plain blue 
curve, the QPI modulation are nicely reproduced with a 
periodicity of ss 0.3 x lQ 1A W/cm 2 as expected. At an 
alignment angle of a = 20o corresponding to the accep- 
tance angle achieved experimentally, the two channel QPI 
pattern, shown in dashed red line, does not differ from 
the single channel one confirming the dominant contribu- 
tion from HOMO (channel X) in the resulting HHG yield. 
However in the case of a — 90o which corresponds to a 
perpendicular alignment of the molecules with respect 
to the laser field, the QPI pattern is clearly shifted by 
7r/2 towards the lower intensity as for the atomic model. 
However, although the periodicity is conserved, the con- 
trast of the QPI is also affected. Therefore we think 
that the shift is a clear signature of the multiple channel 
contribution and that the contrast is a key element to 
identify the laser coupling population transfer between 
HOMO and HOMO-1. 



2. ab-initio TDSE calculation for two electron system 

In order to understand how the laser coupling could 
be disentangled from the multiple channels contribution 
we have performed a more sophisticate calculation based 
on full quantum mechanical calculation on a 2D two- 
electrons model system based on Multi-Configuration 
Time Dependent Hartree method (MCTDHJpSl. The 
method implementation details for the strong field phe- 
nomena can be found irP3 This system corresponds to 
a model molecule with fixed internuclear distance, R = 
3.54 a.u., having a parallel to the molecular axis cou- 
pling between HOMO and HOMO-1 with an energy gap 
of 3 photon energy form the laser fundamental frequency. 
To emphasis the implication of the laser coupling states 
compare to the contribution of multiple channel, we con- 
sidered the case of molecules where nuclear motion is 
slow. Therefore in our calculation, the molecular internu- 
clear distance is frozen. We present two sets of Molecular 
QPI results with this calculation. In one of them only a 
strong 800 nm laser field, polarized perpendicularly to 
the molecular axis, is present. In the other one an ad- 
ditional resonant field, with 3w carrier frequency, polar- 
ized parallel to the molecular axis, is added. In theoreti- 
cal analysis the harmonic emission associated to different 
channels, ionisation from e.g. HOMO and HOMO-1, can 
be resolved^. 

In Fig. 7 we show as an example the results ob- 
tained for the 15th harmonics. QPI pattern for separated 
HOMO (ground state) and HOMO-1 (excited state) with 
and without the resonant field that couples these two 
states. In the case of the absence of the resonant field, the 
HOMO (blue dashed line) and HOMO-1 (green solid line) 
QPI patterns are not shifted with respect to each other, 
since the HOMO-1 state of the ion is populated via the 
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QPI, H15,H 2 (2D), R=3.54 au, 0)^=0.057 au, 10^=30)^, I uv =3.10 13 W/cm 2 




I i I i I i I 
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FIG. 8. QPI pattern for 15th harmonic for a model molecule 
of fixed internuclear distance, R = 3.54 a.u . The lines are 
associated to HOMO without ( blue dashed line) and with 
(black dashed line) the resonant field, and to HOMO-1 with- 
out (green solid line) and with ( red dashed line) the resonant 
field. The strong laser field is polarized perpendicular and the 
resonant laser field parallel to the molecular axis. 

excitation by the recolliding electron onljESl, w hich gives 
no additional accumulated phase. While adding the reso- 
nant field the QPI pattern calculated for HOMO-1 (solid 
red line) presents a shift by quasi 7r with respect to the 
HOMO (black dashed line) due to the additional phase 
caused by the resonant population transfer between the 
cation states. The resonant transition is tuned to ex- 
ceed the transitions via the recollision, by choosing an 
appropriate intensity for the coupling field. Therefore 
the contrast of the HOMO-1 and HOMO populations in 
the presence of resonant field will affect strongly the con- 
trast of the full Molecular QPI signal and is a signature 
of the population transfer. 



IV. CONCLUSION 

In conclusion, we have performed studies that demon- 
strate the potential of Quantum-path Interferences (QPI) 
to probe complex dynamics now in molecules. Our first 
study is experimental and compared short and long tra- 
jectories in randomly aligned molecules (CO2, N2, O2, 
D2, H2). We have clear evidence that the ratio of short 
and long trajectories is sensitive to nuclear motion of 
the cation with attosecond-femtosecond time scale, espe- 
cially in H2 and D2 cases where the comparison is direct. 
We also show that QPI, in randomly aligned molecules 
is governed by the continuum electron phase (as in the 
atomic case) and we highlight the fact that the electronic 
dynamic in the continuum during HHG process is inde- 



pendent on the non-isotropic molecular potential (allow- 
ing the use of strong field approximation for Molecular 
QPI). In a second theoretical study we extend the QPI 
to aligned molecules in order to get inside the electronic 
cation dynamics. We study the cases of multiple channels 
of ionisation contribution as well as population transfer 
due to laser field coupling. In a first simple approach 
based on atomic SFA, we show that the multiple orbitals 
contribution to the Molecular QPI should appear as a 
shift of 7r/2 of the modulation towards low intensities. 
In a quasi-analytical approach based on nitrogen SFA in- 
cluding two ionic states symmetries and their coupling by 
the laser field in perpendicular alignment, the Molecular 
QPI pattern is clearly shifted towards the lower inten- 
sity as for the atomic model but also shows QPI contrast 
losses. This indicates that the shift is a clear signature of 
the multiple channels contribution and the lack of con- 
trast shall allowed to identify the population transfer be- 
tween the ionisation states HOMO and HOMO-1. In a 
third step we calculate the Molecular QPI in a model 
i?2molecule where the population transfer between ionic 
states is allowed by a resonant field. In that case the 
QPI pattern is calculated for each individual channels. 
The HOMO-1 channel reveals a clear shift by quasi 7T 
due to the additional phase caused by the population 
transfer between the cation states. It will therefore af- 
fect strongly the contrast of the full QPI as suggested by 
the quasi-analytic model. 

Hence the QPI method appears to be an efficient 
method to reveal the details of the intra-molecular dy- 
namics in a self-referenced way and applicable to many 
type of molecules. 
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